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Abstract. We present configurations of polarimeters 
which measure the three Stokes parameters /, Q and U 
of the Cosmic Microvifave Background Radiation with a 
nearly diagonal error matrix, independent of the global 
orientation of the polarimeters in the focal plane. These 
configurations also provide the smallest possible error box 
volume. 
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1. Introduction 

This paper originates from preparatory studies for the 
Planck satellite mission. This Cosmic Microwave Back- 
ground (CMB) mapping satellite is designed to be able 
to measure the polarisation of the CMB in several fre- 
quency channels with the sensitivity needed to extract the 
expected cosmological signal. Several authors (see for in- 
stance Rees(1968), Bond and Efstathiou(1987), Melchiorri 
and Vitt orio(|l996| ), Hu and White( |l997| ), Seljak and Zal- 
darriaga( 199^ )), have pointed out that measurements of 
the polarisation of the CMB will help to discriminate 
between cosmological models and to separate the fore- 
grounds. In the theoretical analyses of the polarised power 
spectra, it is in general assumed (explicitly or implicitly) 
that the errors are uncorrelated between the three Stokes 
parameters /, Q and in the reference frame used to 



build the polarised mu ltipoles (Zaldarriaga and Seljak, 
1997] ; |Ng and Liu, 1997| ). However, the errors in the three 



Stokes parameters will in general be correlated, even if 
the noise of the three or more measuring polarimeters are 
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not, unless the layout of the polarimeters is adequately 
chosen. In this paper we construct configurations of the 
relative orientations of the polarimeters, hereafter called 
"Optimised Configurations" (OC), such that, if the noise 
in all polarised bolometers have the same variance and are 
not correlated, the measurement errors in the Stokes pa- 
rameters /, Q and U are independent of the direction of 
the focal plane and decorrelated. Moreover, the volume of 
the error box is minimised. The properties of decorrelation 
and minimum error are maintained when one combines re- 
dundant measurements of the same point of the sky, even 
when the orientation of the focal plane is changed between 
successive measurements. Finally, when combining unpo- 
larised and data from OC's, the resulting errors retain 
their optimised properties. 

In general, the various polarimeters will not have the 
same levels of noise and will be slightly cross-correlated. 
Assuming that these imbalances and cross-correlations are 
small, we show that for OC's the resulting correlations be- 
tween the errors on I, Q and U are also small and easily 
calculated to first order. This remains true when one com- 
bines several measurements on the same point of the sky, 
the correlations get averaged but do not cumulate. 

Finally, we calculate the error matrix between E and 
B multipolar amplitudes and show that it is also simpler 
in OC's 

1.1. General considerations 

In the reference frame where the Stokes parameters /, Q, 
and U are defined, the intensity detected by a polarimeter 
rotated by an angle a with respect to the x axis is: 



la = ^(/ + Q cos2a + [/ sin 2a). 



(1) 



Born and Wolf (1980) page 554, where /, Q, U are noted so, 
si, S2 respectively 



Because polarimeters only measure intensities, angle a 
can be kept between and w. To be able to separate 
the 3 Stokes parameters, at least 3 polarised detectors 
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are needed (or 1 unpolarised and 2 polarised), with angu- 
lar separations different from multiples of 7r/2. If one uses 
n > 3 poiarimeters with orientations ap, 1 < p < n for a 
given line of sight, the Stokes parameters will be estimated 
by minimising the x^: 



- (M - ASf N-^ (M - AS) 



(2) 



/ mi \ 



where M = 



is the vector of measurements, and 



V TO« / 

TV is their n x n noise autocorrelation matrix. The n x 3 
matrix 

/ 1 cos2ai sin2Q!i \ 



A=- 
2 



1 cos 2ap sin 2ap 



\ 1 cos 2a„ sin 2a„ / 



(3) 



relates the results of the n measurements to the vector of 
the Stokes parameters S' — Q in a given reference 

frame, for instance a reference frame fixed with respect to 
the focal instrument. If one looks in the same direction of 
the sky, but with the instrument rotated by an angle ip in 
the focal plane, the matrix A is simply transformed with 
a rotation matrix of angle 21/): 



and the inverse of the covariance matrix of the Stokes 
parameters is given by: 



1 



T 1 

X = A^A = - X 

4 



/ n cos2Q;p J^i sin2ap 

cos2Q;p i(n + cos4ap) i^"sin4ap 



(7) 



\ sin 2ap \ Y.i sin 4Q!p 



cos 4ap) J 



It is shown in the appendix that, if the orientations of the 
poiarimeters are evenly distributed on 180°: 

7r 

ftp = ai + (p — 1) — , p — l...n, with n > 3, (8) 

the matrix V takes the very simple diagonal form: 

/I 0' 

Vo = CT^Xo^S with Xq-^ = - 2 I , (9) 

\0 2 

independent of the orientation of the focal plane, while 
its determinant is minimised. In other words, at the same 
time, the errors on the Stokes parameters get decorrelated, 
their error matrix becomes independent of the orientation 
of the focal plane and the volume of the error ellipsoid 

takes its smallest possible value : ^ ("^) • 

The "Optimised Configurations" (OC) are the sets of 
poiarimeters which satisfy condition (||), (see figure |l]). 
They are hereafter referred to by the subscript as in 



1 

A^AR{ip), with R{ip) = \ cos 2?/' sin2V' | . (4) 

— sin 2-0 cos 2ip 

As is well known, the resulting estimation for the Stokes 
parameters and their variance matrix V are: 



[A^ AT^i A 



and 



V = {A^N-^A 



A^N^'^M, 



(5) 



2. Optimised Configurations 

2.1. The ideal case 

If we assume that the measurements nip (1 < p < n) have 

identical and decorrelated errors 

{Npq =< Snip Sniq >— ctq (5pg), the x^ is simply: 





71/3 



71/4 




71/5 



p=l 



nip — -{I + Q cos 2ap + U sin 2ap) 



(6) Fig. 1. The relative orientations of poiarimeters in "Op- 
timised Configurations" with 3, 4 and 5 detectors. 
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equation (^. The smallest OC involves three polarime- 
ters with relative angle 7r/3. With 4 polarimeters, the an- 
gular separation must be 7r/4, and so on. Note that a con- 
figuration with one unpolarised detector and 2 polarised 
detectors can never measure the Stokes parameters with 
uncorrelated errors, because this would require polarime- 
ters oriented 90° apart from each other, which would not 
allow the breaking of the degeneracy between Q and U. 

If one combines several OC's with several unpolarised 
detectors, all uncorrelated with each other, the resulting 
covariance matrix for the Stokes parameters remains di- 
agonal and independent of the orientation of the various 
OC's. More precisely, when combining the measurements 
of Ut unpolarised detectors (temperature measurements), 
with np polarised detectors arranged in OC's, the covari- 
ance matrix of the Stokes parameters reads: 



V = 



np 



( 















n p 











2 / 



(10) 



where we have introduced inverse average noise levels, gt 
and (Tp, for the unpolarised and polarised detectors re- 
spectively: 



1 



(JT 



, and 



1 



1 



. , . )■ (11) 

^unpolarised / \ *^polariscd / 

Note that the levels of noise can be different from one OC 
to the other and from those of the unpolarised detectors. 

2. 2. A more realistic description of the measurements 

In general one expects that there will be some slight imbal- 
ance and cross-correlation between the noise of the detec- 
tors. The noise matrix of the measurements will in general 
take the form: 

N = a^{l + f3 + f), (12) 
where the imbalance (3 and cross-correlation 7 matrices 

\ 

, Tr(/3) = 0, and 



/3 





1322 
1333 



7 ^ 



V ; : ; 

/ 712 713 

71 2 72 3 

713 72 3 

V : ; ; 



(13) 



are "small", that is will be treated as first order pertur- 
bations in the following, and therefore 

Ar-i = 4(]l-/3-7). 



To this order, the variance matrix of the Stokes parameters 
is easily obtained from Eq. (H): 



y = [x-i + s + g] , 



(14) 



where the matrix X. is given by Eq. (|^) and the first order 
corrections to V, matrices 3 and S, read: 



A X 



(15) 



In an OC, the matrix X^^ takes the simple diagonal form 
Xq^^ of Eq. (^, and the non diagonal parts, "B and S 
remain of order 1 in /3 and 7. For instance, if we consider 
an OC with 3 polarimeters, and polarimeter number 1 is 
oriented in the x direction. 



3 



/ 







/3: 



1 1 

Pii /3ii 

02 2—03 3 03 3—02 



V3 

where /32 2 + /Sa 3 = 1 , and 
„ 4 



02 2—03 3 

03 3—02 2 

V3 
-Pll 



X 



/ 2(71 2-1-71 3-1-72 3) 71 2-1-71 3-272 3 71 2-71 3 

3 3 V3 

71 2+71 3 — 272 3 2(72 3 — 271 2-271 3) 2(71 2— 71 3) 

3 3 73 

71 2— 71 3 2(71 2—71 3) _ 2 '-^2 3 



V3 



V3 



Note that S and S transform under a rotation of the focal 
plane by a rotation R{ip): 



(16) 



Because Vq is invariant, as long as (3 and 7 are small 
the dependence of V on the orientation of the focal plane 
remains weak. 

3. Co-adding measurements 

The planned scanning strategy of Planck goes stepwise: at 
each step the satellite will spin about 100 times around a 
fixed axis, covering the same circular scan, then the spin 
axis of the satellite will be moved by a few arc-minutes, 
and so on. This provides two types of redundancy: every 
pixel along each circle will be scanned about 100 times, 
and some pixels will be seen by several circles, with dif- 
ferent orientations of the focal plane. In this section we 
show, assuming a perfect white noise along each scan, that 
the properties of the error matrix of the Stokes param- 
eters coming from OC's are kept if all data are simply 
co-added at each pixel, whatever the orientations of the 
focal plane. The redundancy provided by intersecting cir- 
cles can be used to remove the stripes induced on maps 
by low-frequency noise in the data streams. An extension 
adapted to polarised measu rements of the method pro- 
posed by Dclabrouillc ( 1998 ) for t he dc -striping of Planck 
maps is studied in Revcnu ct al. ( 199§| ). 
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Here we assume that the noise is not correlated be- 
tween different scans and can thus be described by one 
matrix Ni for each scan, indexed by I, passing through 
the pixel. The is then the sum of the xf over the L 
scans that cross the pixel: 



= J2^Mi - Al SfNr\Mi - Al S). 



(17) 



1=1 



The estimator of the Stokes parameters stemming from 
this x^ is 



S=[J2Af N^' Ai\ ^i' 

/ 1 = 1 

with variance matrix: 

-1 



v=(j:aTn^^a^ 



(18) 



(19) 



In the ideal case, for a given scan, the noise (assumed 
to be white on each scan) has the same variance for all 
bolometers with no correlation between them, although it 
can vary from one scan to the other: 



Ni=aiH, 



(20) 



and one can write the resulting variance combining the L 
scans: 



(21) 



where Xi = Ai Ai, and we have written explicitly the 
rotation matrices which connect the orientation of the fo- 
cal plane along scan / with that along scan 1. Note that 
these matrices are dependent of the position along the 
scan through angle ipi . 

If the observing setup is in an OC, all orientation de- 
pendence drops out and the expression of the covariance 
matrix becomes diagonal as for a single measurement (Eq. 



A 2/100, 
Vo. = i^ 2 U^Xo-^ 
"^ \ 2 ' ^ 



(22) 



where Xq ^ is defined in Eq. (H) and the average noise 
level (TL is defined as: 



err 



(23) 



Of course one recovers the fact that, with L measurements, 
the errors on the Stokes parameters are reduced by a fac- 
tor -\/L. 



More realistically, we expect that the noise matrices will 
take a form similar to Eq. (|l2): 



Ni^ai'{l + (3i+-f,) 



(24) 



If /3; and Ji are small, first order inversion allows to cal- 
culate V (Vl is given by Eq. (|2l|)): 



(25) 



1=1 



If the focal plane is in an OC, this expression simplifies to 



Vl = ^\X^ ■ r ^ 



where 



(26) 



j = R-\'i]h)X„-' Al { j Al Xo-iH(V';)- 



(27) 



The 1/L factor inside the parenthesis in equation ( |26[ ) 
implies that the cross-correlations and the dependence on 
on the orientation t/i; of the focal plan remain weak when 
one cumulates measurements of the same pixel. 

4. The error covariance matrix of the scalar E 
and B parameters 

Scalar polarisation parameters, denoted E and B, have 
been introduced, which do not depend on the reference 



frame ( 


Newman and Penrose, 1966; 


Zaldarriaga and Sel- 


iak, 1997 


). However, the properties of OC's do not propa- 



gate simply to the error matrix of the E and B parameters 
because their definition is non local in terms of the Stokes 
parameters. 

Nevertheless, if the measurements errors are not corre- 
lated between different points of the sky (or if the correla- 
tion has been efficiently suppressed by the data treatment) 
then the elements of the error matrix of the multipolar co- 
efficients aE.im and aB^im can be written in a simple form 
which is given in Appendix B for a general configuration. 

For an OC, the error matrix simplifies further and its 
elements reduce to: 

< SaE i„JaE i,„-^,* >= i (4;—] E^stokcs("fe) 



2 \N., 



X yYlmiflk)* 2Yl'„i'{nk) + -2Yhn{hk)* -2Yl',n'{nk)] 



< 5aE,lm 5aB,Vm* ^ \'^f~ ] E '^Stokos(^fe) 



^ pix 



X [lYlmink)* 2Yl'rn'{nk) - -2Yl„i{nk)* -2Y1' 

where Npi^ is the total number of pixels in the sky, crgtokcs 
is the common r.m.s. error on the Q and U Stokes parame- 
ters, rik is the direction of pixel k and functions ±2^m(^fe) 
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are the spin 2 spherical harmonics. If CTg^Qj^j,;, does not de- 
pend on the direction in the sky, a highly improbable situ- 
ation, then the orthonormality of the spin weighted spher- 
ical harmonics makes the error matrix fully diagonal in the 
limit of a large number of pixels: 



< SaBim Sue I' 



> < SaE,im Sqbm 

> < SaBdm 5aB,v 
An 



1 



<^Stokcs ^11' 



Note that, even for unpolarised data, the error matrix 
between multipolar amplitude is not diagonal unless the 
r.m.s. error is constant over the whole sky (see for instance 
Oh, Spergel and Hinshaw, 1998). 

In the same conditions, the noise matrix of fields E 
and B is also fully diagonal: 



{6E{n) SE{h')) {dE{n) SB{n')) \ _ 



{5B{n) 6E{h')) {SB{h) 6B{h')) 



2 X. 
tokcs "nn' 



5. Conclusion 



In this paper we have shown that, if the noise of the po- 
larimeters have nearly equal levels and are approximately 
uncorrelated, these detectors can be set up in "Optimised 
Configurations" . These configurations are optimised in 
two respects: first the volume of the error ellipsoid is 
minimised, and second the error matrix of the inferred 
Stokes parameters is approximately diagonal in all refer- 
ence frames. This remains true even if one combines in- 
formation from numerous measurements along different 
scanning circles. Such "Optimised Configurations", with 
3 and 4 polarimeters, have been proposed by the Planck 
High Frequency Instrument consortium ( HFI Consortium 
199S) for the three channels where it is intended to mea- 



In such configurations, equation (A.l) becomes: 



i2ai 



i2 n 5a 



i2Sa 



1 



0, 



(A.3) 



The solutions of equation (A.3) under conditions ( |A.2| ) 
reduce to: 

Sa^-, with n > 3. (A.4) 

n 

It is easily seen that conditions ( A.4 ) also automatically 
ensure that 5*4 = and therefore that X becomes diagonal 
and assumes the very simple form: 

Xo = - 1/2 (A.5) 
* \0 1/2/ 

independent of the orientation of the focal plane. Equation 
(^ is the consequence of ( A.5 ) 

The error volume is proportional to the determinant 
of the error matrix V. Therefore, it is minimum when 
the determinant of -X" (equation (^) is maximum. This 
determinant can be written as: 

Det(X) = 

^ [n^ -n\S4\^ -2\S2\\n-\S4\ cos(04 - 2 ^2))] ■ 

(A.6) 

Because the Sk^s are sums of n complex nu mbe rs with 
modulus 1, \Sk\ < n, and it is clear from Eq. (A^) that 

Dct(X) < — , 
^ ^ - 64 ' 

and that the upper bound is reached if and only if 

82^84^ 0. (A.7) 



sure the CMB polarisation. 



Conditions ( [A.2| ) and (A.4) have been shown above to 
imply (A.7), and therefore ensure that the determinant of 
the covariance matrix V is minimum. 



Appendix A: The conditions for an OC 

In this appendix, we show that conditions (||) diagonalise 
the error matrix V of the Stokes parameters and minimise 
its determinant, if the noises in the n polarimeters have 
identical levels and are not correlated. 
We use the notation: 

n 

p=i 

It can be seen from equation (^ that requiring that the 
error on / be decorrelated from the errors on Q and U is 
equivalent to the condition: 



^2-0. 



(A.l) 



This condition can easily be fulfilled in a configuration 
where the angles ap are regularly distributed: 



ap = ai + {p — 1) 6a, p — l...n, 

with n > 3, < Sa < TT, 6a ^ ■n 12 (see text). 



(A.2) 



Appendix B: The general error matrix of the E 
and B multipolar coefficients 

Assuming that the data treatment has removed all corre- 
lation between different directions in the sky, the matrix 
elements of the error matrix of the E and B multipolar 
coefficients are: 

< <Sag,i„ Jag, J ('^') 

pix / 

{Nqq + Nuu){flk)[2Yl^* 2yi'm' + -2yim* -2Yv^>]{nk) 
± {Nqq - Nuu){fik)[2Yl„* -2^1'™' + -2^1™* 2i^'m'](nfc) 
±2iNQu{rik)[2Ylrn* -2Yl'm' - -2Yl„* 2Yv„L']{nk), 



< 6a 



E,lm 



6aB,l'm'* >— -. 



An 



4:\N, 



pix 



{Nqq + Nuu){nk)[2Yl„^* 2Yl>,n' - -2>1™* -251'm'](nfc) 

- i^QQ - Nuu){nk)[2Yl,n* -2Yl'm' - -2Ylm* 2i1'm'](Wfe) 

— 2i NQij{hk)[2Yijn* -2^i'm' + -2^im* 2^('m' ] (l^fe) , 
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where Nifik) is the noise matrix of the Stokes parameters 
Q and U in the direction of pixel k. 
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